1. Introduction 



First order electroweak phase transition has often been discussed [1] within the 
context of the finite temperature effective potential of the order parameter, which 
is the expectation value of the Higgs field, evolving in a thermal background. Such 
phase transitions are highly local phenomena: critical bubbles have a finite size. 
However, effective potentials are averaged over all space and do not describe local 
fluctuations, which may affect the dynamics of nucleation and bubble growth. 
Fluctuations with spatial correlations comparable to the critical bubble size may 
be expected to be important for bubble nucleation and, therefore, we concentrate 
on the thermalization of an on-shell Higgs fleld with zero spatial momentum 
{po = rriH, P = 0). This choice reflects the fact that critical bubbles are typically 
much bigger than the inter particle distance ~ 1/T in plasma. 

Local fluctuations of the order parameter, which one might call subcritical 
fluctuations, are continuously created and destroyed by the self-interactions of 
the thermal background. A natural way to describe them is to consider a 
non-constant order parameter, for which one needs the complete effective action 
for the background field (p{x). The computation of the effective action is in general 
very difficult, but if one focuses on small-amplitude fiuctuations, (j){x) = (j)o + b{x) 
with b{x) small, one may make use of a perturbation expansion in b{x), and to 
order the finite temperature effective action reads simply [2] 

r[</>o,6] = -K//(</>o) + rW(</'o)K^ = o) + ^ I ^b{-k)T(^\<Po;k)bik) , (i) 

where the tadpole F*^^'' = at the high temperature minimum (po = 0, and F''^-* 
is the two-point function for the Higgs field, which can be computed in ordinary 
perturbation theory. 

Two-loop dispersion relations and the evolution of the Fourier modes b{k) 
in Ac^'^-theories has been discussed in detail in [2]. In the present paper we 
investigate how fast the fiuctuations thermalize in the Standard Model above, 
and just prior to, the electroweak phase transition. We assume the standard 
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particle content and one Higgs doublet. The thermalization rate 7 depends on the 
imaginary part of the effective action ImE = — ImF^^^ via 



where a; = a;(fc) is the energy of a given mode. We compute the imaginary part of 
the two-point function at one-loop level for fermion and gauge boson corrections 
(at one loop there is no higgs contribution to ImS in the unbroken phase). The 
imaginary parts are cuts of the relevant one-loop diagrams [3,4], and therefore 
equivalent to processes which may be called decay, inverse decay, absorption, or 
emission of the higgs. For massless fermions and gauge bosons these processes 
are kinematically forbidden at one-loop level. However, the physical external 
states should include thermal corrections which, to leading order in T, are the 
hard thermal loops. They merely change the pole structure of the propagators 
to account for the interactions with the thermal background. Accordingly, there 
appear energy thresholds, which in particular for the fermionic loops are rather 
complicated because of the complicated nature of the dispersion relations. Here, 
the one-loop case is exceptional in that at higher loops there are no thresholds. 

The paper is organized as follows. Section 2 contains a discussion about the 
fermion dispersion relations at finite temperature and different channels, and we 
present a derivation of the fermionic contribution to the thermalization rate. In 
Section 3 we compute the gauge boson loop and estimate also the size of the 
two-loop corrections. In Section 4 we compare the thermalization rate to the 
electroweak phase transition rate, and Section 5 contains our conclusions together 
with a summary of our results. 
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2. Fermions 



2.1 Preliminaries 



The dispersion relation of massless chiral fermions is drastically changed by 
the high temperature correction [5,6]. The leading part of the fermion 
self-energy gives rise to two kinds of excitations of the thermal heat bath called 
particles and holes [7]. Both particles and holes have unusual dispersions relations. 
Particles have a minimum group velocity that is 1/3 of the speed of light, so 
that there are no particle excitations at rest, even though there is a mass gap. 
This phenomenon is due to the collective excitations. Holes have the same 
group velocity for k = but with an opposite sign, and a minimum energy for 
approximately k = 0.408 m/ [8]. All these features combined together make the 
thresholds difficult to determine. The dispersion relations for particles and holes 
are displayed schematically in Fig. 1. 

Before studying the thresholds of the diagram in Fig. 2 we need an expression 
for Im S including all possible channels. We may use the expression in [4] suitably 
modified to account for collective excitations and chirality. One difference between 
the right and left handed fermion propagators is thermally induced masses m^,/? 
which are not equal because the right and left handed fermions are in different 
representations of the Standard Model gauge group. The one-loop resummed 
propagator in the rest frame of the heat bath reads {P^ = (po)P) with p = |p|) 



S{P) = V 



r) + r) 



(3) 



B?{P) 



where 



-(1-75) , ?7 = 2^-^ ' 

(1 + ai(po, V))P^.1^ + ^>l(Po, p)7°, 




(4) 
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and similarly for Iji. The functions aL,R{po-,Y>) ^iid ^l,-r(PO)P) are respectively 
even and odd with respect pq, i.e. they have the properties 

aL,R{-po,v) = aL,R{po,v) , ^L,ii(-po,p) = -hL,R{po,v) ■ (5) 

We have introduced left and right handed fermion masses in Eq. (4). For leptons 
one finds 

^Ir 

where mi and mii are the T = masses of the leptons I and V belonging to the 
same SU{2) doublet. For quarks one finds that 



m 



z + 2m^ + mf + mf, 



2p 



2/2 



(6) 



^Ui,R = 



1 

6 



m 



w. 



+ mi . + m 



di , 



8p 



(7) 



2/2 



and 



di,L 



^^di ,R 



1 



2p 



(8) 



where m^. and m^^ are the T = masses of m and d type quarks of z-th generation, 
respectively, mw = 80.0 GeV, mz = 91.2 GeV and mt is the T = mass of the 
top-quark. The higgs vacuum expectation value is / ~ 247 GeV. For the strong 
coupling constant we adopt the value Qs — 1.2. Note that the left handed masses 
of fermions belonging to same SU (2) doublet are equal due to unbroken symmetry. 
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The imaginary part of the self-energy of the diagram in Fig. 2 is 



ImS/(po,p) = - g.^^26> j ^^^f(^o)e(/co -Po)^sin2(/)fesin2(/)fc. 

X L^{K)R^'{K - P)rjLVRS{L\K))5{R\K - P)) , 




(9) 



where 



sin 20fe = 2 



e/3|feo|/2 



(10) 



ef^l'^o] -I- 1 



and gy is the corresponding Yukawa coupUng. The factors r]L and rjn are +1 and 
— 1 for particles and holes, respectively. They arise when the energies in LP'[K) 
and R^{K — P) are continued analytically using the time-ordered prescription 
^0 ^0 + ie.kQ. For a free massive propagator this would only add ie in the 
denominator. In our case we have 



The factor that multiplies ie is positive for particle solutions and negative for hole 
solutions which gives rise to an extra sign for the hole contribution. Below we 
shall see that this extra sign is crucial in order to get a positive damping rate 
in all channels, whereas the naive propagator (L^(i^) + ze)~^ would lead to an 
instability. 



From the analysis in [3] it follows that there are two possible decay channels for 
a scalar interacting with fermions with two different masses, m\ and 777-2. In the 
case of a scalar at rest (po = fnui P = 0) they are: (1) decay if mn > rrii + m2, (2) 
absorption if ruH < \mi — m2|. The inverses of these processes are also possible 
and they are included in the rates we compute below. Here, and in what follows, 
we denote by mn = mn^T) the T-dependent higgs mass. 

In our case the thresholds are not determined by the k = energy but must be 
derived using the complete dispersion relations. A decay or absorption can take 
place when one of the following conditions is fulfilled 



(fco + zefco , k) = (K) + ze /co 



(11) 



2.2 Thresholds 



(1) mn = a;L,i(k) +a;iij(k) , 

(2) niH + UJL,iO<-) = (^RjOi) , 

(3) iriH + t^fl,i(k) = a;i,,i(k) . 



(12) 
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Here k= |k| and for each energy we have put an index =p,h) to distinguish 
between particles and holes. However, not all combinations of particles and holes 
contribute in each channel because the factor L^R^^ can vanish. This turns out to 
be important in order to obtain a positive thermalization rate. 

Let us denote by E the energy of an excitation, including a sign with reference 
to the arrows in Fig. 2. Then El and Efi have opposite signs in the decay process, 
and equal signs in the absorption process. From Eq. (4) we have explicitly 

L^{EL,i,^)R^'{ER^^,\) = 

[£;L,^(l + aL(EL,,,k)) +6L(EL,i,k)][Efl,,(l + aK(£;K,,-,k)) +6ii(£;Kj,k)] (13) 
- k^ (1 + aL{EL^i, k)) (1 + ajt(Efl,,-, k)) . 

This should be compared with the on-shell conditions for positive energy particles 
and holes 

ujp{l + a{ujp, k)) + h{ujp, k) = k(l + a{ujp, k)) , 

(14) 

a;/i(l + a(a;/i, k)) + h{uh, k) = -k(l + a(u;/,, k)) . 
Negative energy solutions, where uj —uj in Eq. (14), are also possible but then 
the on-shell conditions for particles and holes are reversed due to the properties 
of (^L,R and hL,R in Eq. (5). Using Eq. (14) in Eq. (13) we see that when Ei, 
and En have the same sign, only {i = p,j = h) and (i = h,j = p) give a non-zero 
result. When El and En have opposite signs we are left with {i =p,j = p) and 
(i = h, j = h). This selection rule can be understood physically as a consequence of 
the conservation of angular momentum. The chirality-helicity relation is reversed 
for holes [7], and because the fermions in scattering processes under consideration 
always have different chiralities, absorption is possible only if they have same 
helicities, and direct decay is possible only if they have opposite helicities. 

For a Higgs particle at rest the angular integrals in Eq. (9) become trivial and 
the 5-functions can be used to perform the ko and k integrations. One should note 
that the (5-function may be rewritten as 

5{L^) = J2 WL,i{k)[5{ko-u;LM) + 5{ko + u;L,i{^))] > (15) 

i=p,h 
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where the energy factor WL,i is different from the usual l/2a;(k). In fact, we have 
instead 



dko 



(16) 



ko=U>L,i 

and similarly for Wr. Combining the expression in Eq. (13) for L^R^^ with the 
factors Wl and we find the rather simple expression 

2«,-k2)(c.|,^.-k2) 



L^RI'WlWr = 



(17) 



which will be useful later. We note that L^^W^WlWh is always negative. 



2.3 Decay rates 



Let us start with the absorption channel, which opens up when ttih is small 
enough. Its existence is determined by the last two conditions in Eq. (12). In 
this case El and Er have the same sign so that only p h and h ^ p processes 
contribute. By virtue of the dispersion relations (see Fig. 1) there is a channel for 
small niH even if m_L = m_R, because when coh < (^p the higgs may be absorbed by 
a hole in the heat bath to produce a particle. Using Eq. (17) the imaginary part 
for the process H + fu^h — > /^^p can then be found to be 

Im Lp) = ^Vl,pVrA' ^""''^ 'a^IaIi " [^^-^(^^ - ^^'^(^)] ' ^^^^ 

where k is determined by the relation 

rriH + (^R,h{^) = ^L,p(k) , (19) 

and 

1 , , 

= e/3u^.,H + 1 ' (20) 

are the fermion distribution functions. We find that ImS in Eq. (18) is negative, 
as it should, leading to a positive decay rate, but the reason is different from the 

case of ordinary massive fermions [3]. Here, the matrix element L^R^ is negative 
in all channels, in contrast to the case in [3] where the signs are different for decay 



and absorption. On the other hand, we have an extra factor of r]L,ir]Rj which 
then helps to give the correct sign. We note that this works because only p ^ h 
and h ^ p processes contribute in this channel. To obtain the complete fermionic 
Im S for absorptions, the contribution from the process -ff + /l,^ — > fR,p should be 
added to Eq. (18) so that Im T,abs = Im T,{Rh, Lp) + Im T,{Lh, Rp). In the minimal 
Standard Model the difference between rriL and mR is so small compared to the 
particle-hole difference that we could safely approximate rriL — ruR, which would 
simplify the expression for ImS^bs- (See, however, the discussion below when the 
mass difference is not neglected.) 

From Fig. 1 we see that the threshold is at a finite ki and two solutions around 
k = ki exist for niR slightly below the threshold. When niH decreases one solution 
goes to zero, if we neglect the mass difference, and therefore its contribution to 
Im S goes to zero. The other solution goes to infinity but the residue of the hole 
propagator, which is proportional to a;^ — k^, goes exponentially to zero and thus 
Im Safes — > when mn — 0. 

If we do not neglect the mass difference in the absorption channel the threshold 
structure becomes more complicated. Let us assume that > rriR. Then 
in the H + /r^^ f^^p channel (Fig. la) one of the solutions disappears when 
rriH < tul — rriR. In the H + f^^^ — > fR p channel (Fig. lb) the dispersion curves 
cross at some point k2. The threshold is determined by the largest diff'erence to 
the right of the crossing point (ki in Fig. lb). As uih decreases there appear 
two solutions to the right. In the limit uih — ^ 0, the large k solution disappears 
since the hole residue goes exponentially to zero. We are then left with a solution 
which for a small mass difference is located at k2 = |(mi, — itir). When uih —>■ 
the thermalization rate can be found from Eq. (18) to be 

„, /rp ^ ^9y ( i^R + V - mL y 1 

IGtt V ^rriRmL ) \ T ) 1 + cosh[/3(mK + m^)] ' ^ ' 

There is also a third absorption channel, H + /^^ p fL,hi which opens for 
niH < TTiL —iTtR and k is located between and k2. The thermalization rate in 
the limit mn — >^ is again given by Eq. (21). 



The direct decay channel is simpler. The signs of El and En are opposite, and 
the higgs decays either into a pair of particles or into a pair of holes. Therefore, 
the product tjltjr is +1, leading again to a positive decay rate. Summing over 
both decay channels we find 



ImSdec= >^ k , 2 A 2 [i- - nL,i - riR^i) , 

TT Ami Am- 



(22) 



i=p,h 

and k is in this case determined by 



rriH = WL,i(k) + ujR,i{k) . (23) 

The ordinary T = contribution is included in this expression. The threshold 
for holes is lower than for particles since the minimum energy for holes 
is about 0.928 m/ [8]. There are also two solutions in this channel when 
1.865 m/ < niH < 2 m/. As tuh increases the corresponding k increases and the 
hole residue goes exponentially to zero so that the particle channel dominates as 
soon as ruH is above the particle threshold. 

The full fermionic contribution is thus finally given by Im S/ = Im Sa;,s+Im Ti^ec- 
In practise the direct decay channel is in many cases irrelevant, because its 
existence would imply that niHiT) > T. At least in perturbation theory m^ is 
expected to be less than T in the high temperature limit. We shall return to 
Im S/ in Section 4, where we also plot it for the top quark. 

3. Gauge bosons 



3.1 One- loop case 



At T = the imaginary part of the diagram of Fig. 3 is zero for an on-shell 
massive Higgs particle and a massless gauge boson for kinematical reasons. At 
finite temperature thermal masses are induced for both the Higgs and the gauge 
boson so that an absorption channel opens when 2m//(r) < mw{T). In particular, 
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this condition may be fulfilled close to the critical temperature of a second order, 
or a weakly first order, phase transition where the higgs mass becomes small. 

The self-energy of a gauge boson is a rather complicated function of the 
momentum [5,6]. We may, however, still use the standard expression for the 
imaginary part of the self-energy, which can be found in e.g. [4], if we take into 
account the fact that the higgs couples to gauge bosons through a derivative 
coupling, and replace the 5-functions with the ones containing the thermal on-shell 
condition. With the high T expression for the gauge boson propagator* [5,6] we 
get (K^ = (A;o,k),k=|k|) 

ImEgb(po,p) = 

f d^K 1 



sinh2^, / -h2^,sinh2V. ^((p - k)^ - m%{T)) ^^^^ 



where 



and 



e/?IPo|/2 

sinh2^„ = 2^ . (25) 



.2 ■ 4 



for SU{2) , 



= <; (26) 
^'2 for U{1) . 

The gauge boson self-energy has a transverse and a longitudinal part, which are 
given by 



7rT{ko,k) = 



7rL{ko,k) = 



{nw + 1) 



k2 ^ k2 ^ 2k 



3 k2 ^ 



1-^ln 
2k 



kr 



ko-k 



(27) 



These self-energies include the contributions from gauge boson, fermion and scalar 
loops. We use the notation nw = {Nf+AN)/A, where ATj = 2 x (# o/ generations) 

* We are using the Feynman gauge throughout the paper. 
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and N = 2 for SU{2) gauge fields. For U{1) gauge field g is replaced by g' and 
nw by = ||A^/. We define also the longitudinal and transverse gauge boson 
masses by M£ = 77^(0;, k = 0) and M|, = TTTioo, k = 0), respectively. 

At one loop the high T expression for the higgs mass in the unbroken phase is 
given by 

2 (rj.. rnljOi) 6mit. + 3m| + 6m| + m|(0) 3 

mH{^) = 2 ^ 12/2 ^ ' 

where 777,^(0) is the zero temperature higgs mass. The longitudinal masses of 
SU{2) and U{1) gauge bosons are given by 



while the transverse masses are 



^2 _ nr + 1 4(m| - m^) ^2 



(29) 



9 r 



2 _ -2 4m^ ^2 
2 _ ny -2 4(m| -m^) ^2 



(30) 



It is the scalar loop that gives rise to the difference between the longitudinal and 
transverse masses. 

For a Higgs particle at rest {po — mH,P — 0) only the longitudinal part 
contributes and the angular integral becomes again trivial. We obtain, as in 
Section 2, an energy factor from the residue at the pole, given by 



d 



^ Ml -7rLMk),k) _ ^^^^ 



dko 

The absorption channel opens up when the higgs mass is so small that it can 
annihilate with a higgs from the heat bath to produce a gauge boson. In that 
channel Im S is given by 

IniS,5 = — ^ffi==ml,(T)k^[nH(k) - nz.(k)] , (32) 
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where k is determined by solving the on-sheU condition of the gauge boson relation 

a;i(k) = k2+7rL(a;^(k),k) , (33) 



and the energy conservation condition 



mniT) + ^Jk^ + m%{T) = ^k^ + 7ri,(a;z.(k), k) . 



(34) 



The thermal distribution functions are 



n{k)H,L 



1 



(35) 



with 



^k2 + m2,(T) , u;l = V^ 



(36) 



Equation (32) expresses the imaginary part for any component of the Higgs 
doublet; they are all equal because of the unbroken gauge symmetry. 

We are mainly interested in the thermalization rate close to the phase transition 
where mniT) is expected to be small. In the limit mn — > the thermalization 
rate due to SU{2) gauge bosons is given by 



This result can be found by using the asymptotic formula [5] ttl = 
4k2exp (-6k2/(nvi/ + l)fir2T2) which is valid for 3k/y/nw + 1 fifT > 1, together 
with Eqs. (33) and (34). One solves k iteratively as a function of mn, whence one 
obtains Eq. (37) from Eq. (32), and the expansion is valid when mn ^ 10~^T, 
corresponding to ^gb/T < 1.9 x 10~^. 




{nw + l)g^ , T 



3/2 



6 niH 



(37) 




3.2 Two- loop case 



At two-loop level the variety of relevant processes is much larger than at 
one loop, leading to tens of different diagrams. We now estimate the order of 
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magnitude of two-loop processes using the two representative graphs given in Fig. 
4. The purely scalar "rising sun" diagram of Fig. 4a has been calculated in [2]. 
The imaginary part of the two-point function for an on-shell complex higgs with 
a vanishing external three momentum is, in the leading high T approximation, 

ImSn = . (38) 

2567r ^ ^ 

We may use this result to estimate the imaginary part of Fig. 4b which then reads 

ImS = -T^T^ (39) 
2567r ^ ^ 

for SU{2) gauge bosons. The U{1) gauge boson has a similar expression but 
with replaced by 4(5r')^. In the derivation of Eq. (39) the mass difference 
of gauge bosons and Higgs bosons was neglected; formally, the omitted terms 
are next-to-leading corrections in terms of the gauge couplings. The two-loop 
thermalization rate is thus 

^if/T^^^. (40) 

When niH decreases increases. Therefore, for very small higgs masses the 
two-loop processes could actually dominate over the one-loop ones, for which the 
rate goes to zero as tjih — > as is evident from Eq. (37). 

There are, however, certain issues which make our two-loop estimate Eq. (40) 
uncertain. First, we expect ImE(p^ = 0) to be zero (which it formally is in the 
full two-loop expression in [2]) but infrared dominance leads to the constant in 
Eq. (38), which then gives a diverging 7^ when ttih — 0. Secondly, the dispersion 
relation with a complex E reads 

2 2,2 I / 
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u> 

which modifies the simple relation o;^ = p^ +m|f(T) when m|^(T) < |ImE|. (Note 
that ReS has been included in the definition of the renormalized mass tuh-) 
Also, when propagators have complex poles the 5-functions in Eqs. (9) and (24) 
should be replaced by Lorentz functions. Taken all together, we conclude that 
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the simple comparison with the "rising sun" diagram gives a indication that the 
two-loop contributions may be comparable to the one-loop ones, although the 
actual magnitude remains somewhat uncertain. 



4. Thermalization rate vs. phase transition time 



As an application, let us now turn to consider a first-order electroweak phase 
transition. There has been a number of attempts [1] to derive the form of 
the appropriate effective potential for the Higgs field in perturbation theory, but 
lattice studies [9] seem to indicate that there also exists a large non-perturbative 
component. Therefore, to be definite, we shall merely assume that the free energy 
of the neutral Higgs field for ~ and near the critical temperature Tg can be 
expressed through the expansion 

V{cj>) = Ik{T' - T^)ct>' - laT4>^ + \xcj>' + .... (42) 

where k, ck, A and Tq are (unknown) parameters with = To/(l — 2q;^/9Ak). 
In principle, they need not be equal to the corresponding parameters at T ~ 0. 
Thermodynamical properties of Eq. (42), as well as bubble nucleation, has been 

discussed in [10]. Nucleation starts at Tj where Tq < T/ < Tc, and the time scale 
after which the major part of a given volume is filled by bubbles of new phase can 
be estimated to be 

1 ^A^ 

^growth = cvTTT ~ 3~ 7~i 77 ' y'*'^) 
S'{tf) In^ (M^/T4) 

where S{t) is the action, and in the small supercooling limit A = 
0.0144 Q!^A~'''/^K~^(Tc/To)'^. Here tc is the time corresponding to Tc and is related 
to the Hubble rate hy H = {Sir^g^/QOy/'^T^ /Mp = l/2tc. The size of the critical 
bubble, i?c, depends on the parameters of the potential Eq. (42) but is typically 
larger than the correlation length [11]. 

Thermal fluctuations generate configurations which may locally affect the 
dynamics of the phase transition, and for that purpose it is important to find 
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out the lifetime of these fluctuations. In particular, fluctuations larger than the 
correlation length, if sufficiently stable, could be important for bubble formation 
and growth. 

In the previous Sections we have computed the thermalization rate for 
perturbations with small amplitude {b{x)'^ < {p"^ + rn'jj) / X) and large spatial size 
(p ~ 0); these correspond then to large subcritical "bubbles". Since the precise 
T-dependence of mn close to Tc is not reliably described by the loop expansion, 
we prefer to parametrize 7 as a function of mn rather than T. When varying 
uih/T we can consider T to be fixed and vary the zero temperature higgs mass 
mniO), which is not experimentally known. If we rather want 7 as a function of 
T, we would have to take into account the T-dependence of rriH, and the fact 
that the vertical scale in Fig. 5 depends on T. 

The largest fermionic contribution to the imaginary part comes from the 
top quark intermediate states, and in Fig. 5 we have drawn the fermionic 
thermalization rate, assuming rrit = 135 GeV and mi, = rriji, together with the 
gauge boson contributions. The thresholds are such that the one-loop rate is 
zero when 0.25T < m//(T) < 1.15T. (The treshold structure for the bottom quark 
is different from that of the top quark and would actually yield a non-zero 
thermalization rate in parts of the region 0.25T < mniT) < 1.15T, but with a 
relatively small rate.) In contrast, the two-loop rate does not have any thresholds 
so that thermalization will in fact take place for all values of mn- Comparing 
with Fig. 5, we find that the two-loop estimate Eq. (40) would become bigger 
than the one-loop rate when uih/T < 5 x 10~^. 

From Fig. 5 we see also that the thermalization rate is roughly 7 ~ 10~^T 3> H 
in the region where uih is small. Thus compared with the Hubble time, 
thermalization is very fast. However, the phase transition may also be very fast. 
This is so especially when the amount of supercooling is small. In that case 
typically mn <S T during the transition and we may safely approximate tf ~ tc- 
Comparing thermalization rate with the growth time we find that 




(44) 
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At Tc we may take for reference values the one-loop results k = 0.23 (with 
rrit = 135 GeV) and a = 0.028, and for A an arbitrary value 0.006 to satisfy the 
baryon-number retention bound [12] (in A0^-theories there are actually indications 
[13] that A could be driven to very small values at Tc). For these values ^ = 0.35. 
With Tc ~ 100 GeV, we then find that 7 tgrowth — 5 x 10^. In order that 
7 ^growth ~ 1 we would need A ~ lO"^'^, or extremely small a and phase transition 
which is almost second order. Although in principle this might be possible, we 
conclude that it is unlikely that spatially large fluctuations can remain stable 
during the electroweak phase transition. 



5. Summary and conclusions 



We have calculated the one-loop thermalization rate of the Higgs field in the 
high temperature phase in the Standard Model. A fiuctuation of the Higgs field 
thermalizes through scattering processes with itself, with the gauge bosons and 
with the fermions. It is necessary to resum the hard thermal loops for the gauge 
boson and the fermion self-energies in order to get any kinematically allowed 
processes at the one-loop level. This changes the dispersion relations and induces 
thermal masses in such a way that the Higgs particle can annihilate with a higgs 
from the heat bath to create a gauge boson, be absorbed by a fermion, or decay 
into a pair of fermions. The inverses of these processes are also possible. By using 
resummed propagators for gauge bosons and fermions we consider the external 
states to be the effective particle excitations in the heat bath, including collective 
excitations. 

It is the existence of the two fermionic branches at high T, particles and holes 
with opposite chirality-helicity relations, that allow for a Higgs particle to be 
absorbed by fermions in the heat bath. The unusual dispersion relation for holes 
makes the determination of energy thresholds complicated and it has to be done 
numerically. To get a positive thermalization rate one has to use propagators with 
the correct analytical continuation, ko ^ ko + ieko, and not just add ie in the 
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denominator. Only processes which conserve hehcity are possible, and this is also 
essential in order to get a positive thermalization rate. 

Two Higgs particles can annihilate and produce a gauge boson when the higgs 
mass is small enough. There are, of course, also higher order processes that do not 
have any thresholds. In order to make a crude estimate of the two-loop processes 
we estimate them by a two-loop gauge boson "rising sun" diagram. We found that 
for the Standard Model parameters there is a possibility that the two-loop effects 
are of the same order of magnitude as the one-loop ones. 

The thermalization rate is compared with other time scales in Section 4. 
Two such scales of interest are the Hubble expansion and the electroweak phase 
transition rates. The time it takes the bubbles of new phase to fill up the universe 
is to a large extent model dependent but can be estimated in a phenomenological 
model [10]. We found that in general the thermalization rate is several orders of 
magnitude faster than both the Hubble rate and the phase transition rate. This 
does not mean that thermal fiuctuations are unimportant for the dynamics of 
the phase transition. Rather, it justifies time averaging over fluctuations when 
comparing them with the phase transition rate. 
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Figure Captions 



Fig. 1. Schematic finite temperature dispersion relations for fermions: (a) left 
handed particle and right handed hole; ki marks the maximum difference between 
the particle and hole energies; (b) left handed hole and right handed particle; k2 
denotes the crossing point where the energies are equal. 

Fig. 2. The self-energy diagram for the fermionic thermalization rate. The arrows 
indicate the flow of four-momenta. 

Fig. 3. The self-energy diagram for the gauge bosonic thermalization rate. The 
arrows indicate the flow of four-momenta. 

Fig. 4. Two two-loop diagrams contributing to the thermalization rate: (a) pure 
scalar "rising sun"; (b) "rising sun" with gauge bosons. 

Fig. 5. One-loop thermalization rates: SU(2) gauge boson contribution (dotted 
line), U(l) gauge boson contribution (dashed line) and top quark contribution for 
the absorption channel (a) and the decay channel (b) (solid lines). 
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Abstract 

The thermalization rate for long wavelength fluctuations in the Higgs field is 
calculated from the imaginary part of the finite temperature effective action in the 
unbroken phase of the Standard Model. We use improved propagators including 
a resummation of hard thermal loops. The thermalization rate is computed at 
one loop level, but an estimate of the two-loop contribution appears to give an 
indication that they are comparable to the one-loop result for small thermal higgs 
mass. We show also that the Higgs field fiuctuations are likely to thermalize very 
fast compared with the electroweak phase transition time. 
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